Nernst's heat theorem for Casimir-Lifshitz free energy 
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By regarding the Lifshitz expression for the Casimir free energy on the real frequency axis rather 
than the imaginary Matsubara frequencies as is customary, new light is shed on the ongoing debate 
regarding the thermodynamical consistency of this theory in combination with common permittivity 
models. It is argued that when permittivity is temperature independent over a temperature interval 
including zero temperature, a cavity made of causal material with continuous dispersion properties 
separated by vacuum cannot violate Nernst's theorem (the third law of thermodynamics). The 
purported violation of this theorem pertains to divergencies in the double limit in which frequency 
and temperature vanish simultaneously. While any model should abide by the laws of thermody- 
namics within its range of applicability, we emphasise that the Nernst heat theorem is a relevant 
criterion for choosing amongst candidate theories only when these theories are fully applicable at 
zero temperature and frequency. 
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Since Bostrom and Sernelius first predicted the exis- 
tence ofjarge thermal corrections to the Casimir force 
in 2000'1,], controversies over the thermal behaviour of 
this effect, which in its most typical embodiment may be 
seen as the attraction between macroscopic objects due 
to zero-point fluctuations of the quantum vacuum, have 
been extensively covered in the published hterature. The 
use of the Drude model to describe dielectric permittivity, 
employed in [l| was soon criticised on thermodynamical 
grounds Q . The reason was that in the case of a perfect 
crystal lattice, when all dissipation is due to scattering 
of electrons on thermal phonons, the Casimir free en- 
ergy as calculated with the Lifshitz formula appears to 
violate the Nernst's heat theorem which states that en- 
tropy should vanish as T — » 0. The Drude model was 
defended by other authors [1, |3) Hi who argued that since 
the Drude model offers better description of impure met- 
als, and since real metal samples always have impurities, 
the Drude model must be employed. It was shown in 2003 
and recently in a more extensive treatment that the 
free energy with Drude permittivity is quadratic in T for 
small temperatures when impurities are present. No con- 
sensus has yet been reached on the important physical 
question of why Casimir force predictions for the perfect 
lattice model, important in solid state physics, differs sig- 
nificantly from those pertaining to real metals with a very 
small but nonzero concentration of imperfections. 

Recently, a somewhat analoguous problem was brought 
forth for dielectrics with a^small conductivity for finite T 
which vanishes at T = 0[8|. While the purported viola- 
tion in the case of Drude metals referred to the transverse 
electric (TE) mode, this time the bother appears to be a 
discontinuity in the transverse magnetic (TM) Fresnel re- 
flection coefficient giving rise to nonzero entropy at zero 
temperature. The problem was recently argued to extend 
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to insulators, intrinsic and lightly doped semiconductors 
as well as Mott-Hubbard semiconductors and indeed the 
permittivity contribution from Debye rotation of molec- 
ular dipoles[§|. 

According to the argument presented herein both the 
Drude model and the dielectric permittivity model with 
conductivity included belong to a group of permittivities 
which cannot violate Nernst's theorem when permittivity 
can be regarded as temperature invariant in a range of 
temperatures near and including T — 0. 

While the present paper does not aspire to solve the 
physical question of how to take into account the pres- 
ence of a small conductivity in dielectrics when substi- 
tuted into the Lifshitz formula, it seeks to illuminate the 
ever recurring question of thermodynamical consistency. 
It has previously been shown [T^ by use of the Euler- 
Maclaurin (or equivalently Abel-Plana) formula that for 
Fresnel reflection coefficients which are continuous func- 
tions of imaginary frequencies in the limit T — > 0, 
Nernst's theorem is satisfied. An exploration of Casimir 
entropy in the formalism of surface modes was under- 
taken independently of this work by Intravaia and Henkel 
[Tl| whose conclusions accord with ours. By a method of 
summation of the eigenmodes of the vacuum between two 
plates they demonstrate that Nernst's theorem can be 
broken between metal plates only for temperature depen- 
dent relaxation such as in a perfect and infinitely large 
metal lattice. 

This paper demonstrates a similar result using the 
real frequency Lifshitz formalism between plates of a 
generic nonmagnetic materials whose permittivity satis- 
fies a small set of criteria. The real frequency formalism 
is more complicated and less elegant, but with the ad- 
vantage of a more direct physical interpretation. Finally, 
a discussion of the physical implications of the mathe- 
matical limits involved is given. In particular we empha- 
sise the importance of assessing when Nernst's theorem, 
which concerns zero temperature, can be used to inform 
finite temperature physics. 
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I. INTRODUCTION: FREE ENERGY AND 
ENTROPY AT REAL FREQUENCIES 

The Lifshitz expression for the free energy per unit 
transverse area between two identical dielectric plates 
separated by vacuum is in general of the form 

^(a,r) = ^ dwcoth (^2^^ 5m{<?!>(tj,T)}, (1) 

where lot = k^T/h and (l){u!,T) is the zero temperature 
integrand 

^ poo TM 

^{uj,T) = ^ / dk^k^Y. ^riD,{u;,ck^,T) (2) 
wherein 



q=TE 



fTE 



q = 1 - exp(-2Koa); (3a) 



Kq — K 
K0 + 



rTM 



eKo — K 



(3b) 



We have assumed fJ- — 1 everywhere for simplicity, e is 
the permittivity of the dielectric relative to vacuum and 

KQ = (fci - c^Vc2)i/2; n^ikl- ec^Vc')'/'. (4) 

We assume that e does not depend on transverse momen- 
tum, thus neglecting any nonlocal effects and furthermore 
that it is a generalised susceptibility and obeys causality, 
which implies in particular that [l3| 

1. e(— w*) = e*(w) 

2. |3m{e(ci;)}| > on the entire real frequency axis 
except at = where it may be undefined. 

In general e is also temperature dependent, making for 
the temperature dependence of 0(cli,T). The complex 
conjugate is denoted with an asterisk and we will consider 
only real frequencies henceforth. One might furthermore 
impose the physically reasonable demand that 

3. e(w) is continuous and |e(ti^)| < oo for all real fre- 
quencies except possibly a; = 0. 

The function obeys the symmetry property (j){—uj) — 
(j)*{co) for real frequencies [11], hence the real part of (/)(aj) 
is even with respect to oj whilst its imaginary part is odd. 
This allows us to write ^ in a form which makes the 
mathematical discussion in the following somewhat more 
transparent. Since both 5m(/) and coth(ci;/2a;7') are odd 
functions of u, the integrand of ([1]) is even and we can 
let the Lo integral run from — oo to cxd and divide by 2. 
Adding the real part of (f) by substituting Smt/) (jy/i will 
make no difference, since it makes for an odd integrand 
term which vanishes under symmetrical integration, so 



(5) 



is equivalent to H]). 

Assume for the moment that e, and hence (p, is invari- 
ant with temperature over at least a finite range of small 
temperatures including T — Q. In this case the temper- 
ature dependence of J^(T) can be treated very simply 
when T is in this range, since the T dependence now sits 
only in the factor coth(cj/2(jjT). 

From thermodynamics the Casimir entropy in the cav- 
ity, S, is given as 



d_ 

' dT 



(6) 



so if one were able to interchange integration and differ- 
entiation with respect to T, one could write 



5*= - 



l|j.^(.)Aeoth(£-) 

uj exp{uj / lot) 



iksT^ J- 



dL0(j}{L0)- 



eyi\){uj / lot) — 1]^ 



(7) 



For any finite lo the integrand of ([7]) vanishes extremely 
fast, as exp(— /i|cl'|//cbT'), when T ^ 0. This demon- 
strates the finding of Torgerson and Lamoreaux jl4i] that 
temperature corrections are important only for frequen- 
cies below LOT^ which is a very low frequency even at room 
temperature (~ lO^^rad/s). At a; = (and finite T) the 
rightmost fraction in ([7]) has a simple pole, yet only the 
imaginary part of 4'{uj) contributes to ([7|), which is zero 
here since 3m{(/)} is an odd function of w, removing this 
pole. In the limit T — > 0, thus, entropy vanishes as it 
should and the third law of thermodynamics is obeyed 
(there are subtleties pertaining to the TM mode as will 
be discussed in the following). 

Two questions arise from this consideration: Under 
which circumstances may differentiation be interchanged 
with integration? And what happens if one or more pa- 
rameters of t{Lo) are temperature dependent all the way 
down to zero temperature? 



II. TEMPERATURE INDEPENDENT e{Lo) 

We will treat the first question first. Leibnitz' integral 
rule for improper integrals. 



d f°° f°° d 

— I dxf{x,y) = / dx— 



lix.y), (8) 



is always valid when ([T^ §4.44) 

• f{x,y) and df{x, y)/dy are both continuous on x S 
[xq, 00) and the relevant interval of values of y, 

• the integral on the left exists, and 

• the integral on the right converges uniformly. 

The generalisation to integrals with both limits infinite 
is trivial. 
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To make our considerations more concrete, let us con- 
centrate on some permittivity models which are in com- 
mon use: 

= 1 - (9) 



Lu{uj + iv) 



e{uj) = 1 + 



1 



- (10) 



of which the former is the Drude model for metals, and 
the latter describes a semiconductor. Here ojp is the 
plasma frequency, v the relaxation frequency, eq the vac- 
uum permittivity and tooi 7 and ujq material parameters. 
a is the DC conductivity of the semiconductor. Some 
common models which obey [T] but not [2] are 

6H = 1 - (11) 



6(C.) 



1 



1 



1 - wV^o^ 



(12) 



the plasma model[3a| for metals and a model of dielectrics 
with 5 function dissipation ai lo — loq [17| . Notice that 
dSl) and ^ obey criteria [DEI 



A. Propagating and evanescent waves 

One notices that we seem to run into trouble with the 
continuity criterion at \uj\ = cfcj^, where kq = 0, since, 
when regarded as a double integral, ([5]) seems at first 
glance to imply integration across the lines lu = ±cfcj^, 
which would cause trouble with continuity: one sees from 
((3a| that Dq = for kq = 0, hence the real part of 
LnZJq is undefined and the imaginary part turns out to 
be discontinuous as these lines are crossed. 

The problem can be avoided, however. Let us define 
(3 = k^c/uj for short. For positive frequencies, /3 = 1 is 
the limit in which the electromagnetic fields in the cav- 
ity travel parallel to the plates and become evanescent in 
vacuum as the P — 1 barrier is crossed, a limit whose dis- 
continuous properties are physically obvious: the waves 
just on the propagating side {ck± just smaller than uj) 
travel through the system just gracing the surfaces, while 
the fields on the evanescent side stay on the surfaces; they 
are qualitatively different phenomena and the transition 
from one to the other can be expected to be discontin- 
uous |37[. Negative frequencies have no direct physical 
meaning, hence the terms "propagating" and "evanes- 



cent" must be understood in a mathematical sense here, 
defined by |/?| < 1 and |/3| > 1 respectively. 

In the original Lifshitz paper [l3| , the integral is split 
automatically into propagating and evanescent parts by 
substituting p = iKQc/uj. Propagating contributions cor- 
respond to integrating p from 1 to and evanescent to 
an integral from iO to zoo, thus avoiding the problem. 
We notice furthermore that the issues related to \f3\ — 1 
occur for any choice of e{w) hence can have nothing to 
do with the problems with Nernst's theorem which all 
concern particular permittivity models. 



B. Continuity 

In the classical treatment by Casimir, the vacuum en- 
ergy shift was found by summing over the cavity modes of 
the system[l^, a method developed further by van Kam- 
pen et alia by use of the so-called generalised argument 
principle ^1^ and elaborated by Barash and Ginzburg[20j. 
The normal modes of the cavity solve the characteristic 
equation of the set of electromagnetic boundary condi- 
tions which reduce to the equation Dq — 0. At these 
frequencies >p{uj) would have poles which would cause 
trouble with continuity. 

With permittivity models such as (O and (fTU)) where 
dissipation is included (i.e. e{uj) has a nonzero imaginary 
part), D{uj,k±) = has no real- frequency solutions [3g| 
except possibly lu — since rq ^ 1 everywhere. The same 
is the case with (1121) if an imaginary term is inserted in 
the denominator as in (|10p (otherwise \rq \ = 1 at co — luq) . 
The real-valued permittivities as given in ([TT|) and (fT^ 
cause TTM to diverge where eno + k = 0, however, in 
transgression of the continuity criterion. 

We conclude that the continuity of (j>{uj) is ensured for 
all a; ^ so long as is finite, continuous and / 1, 
sufficient criteria for which are that e(cL') satisfies criteria 

m-m 

What remains is the point a; = 0. A priori, this is the 
interesting limit, since when T is very small, the coth 
function in ^ differs from unity only very close to zero 
frequency. As is well known, reflection coefficients are 
occasionally ill defined in the limit where lu and ck± both 
approach zero, as is the case for the Drude model TE 
reflection coefficient, for example, is always bounded, 
however, so the integrand of 4'{uj) approaches zero in this 
limit due to the factor k± stemming from the isotropic 
infinitesimal d'^k± = 2Trk±dk± for any /3 7^ 1. Hence 
(p{u>) is continuous for all uj if e{uj) obeys criterion [31 

A more serious problem is caused by the simple poles 
of coth(a;/2wT) and its T-derivative at w = 0. As argued 
previously, the imaginary part oi(f>(to) is zero at = 0, so 
the integrand of ([7|) does not diverge, but is in some cases 
finite in this limit. For sufficiently small co and finite a, 
(HI) and ((lOl) both have the form e A + iB/to where A 
and B are constants, while if cr = (fTU)) instead has the 
form e A + iBlu. In both cases the imaginary part of 
r^g falls off quickly, as cu^ and cu^ respectively, but when 
A 1, 3m{r^,^} decreases only linearly. One easily 
verifies that with respect to w, 3m{0(a;)} tx 3m{r^(a;)} 
to leading order, hence the TM mode term of 4>(lo) is 
proportional to co in the above mentioned cases. 

To see how this is troublesome, consider the functions 



E/2t _|_ g~x/2t 



I{x,t) 



e 

^2 1 1 



t2[e^/t _ 1)2' 



(13) 
(14) 



which are essentially (up to a constant factor) the inte- 
grands of ([5]) and ([7]) respectively when 0(w) oc w; here 
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FIG. 1: Plot of ^ and ^ as functions of x and t. The 
functions H and / are essentially the integrands of ([SJ and 
^ when </>(cj) oc uj. 



X and t are suitably non-dimensionalised frequency and 
temperature respectively so that uj/ldt — x/t. Notice 
that dH/dt = 21. Equations (HH) and HI]) are plotted 
in figure [1] for x and t. If we define H{Q, 0) to be its 
limiting value 0, the integrand of ([5]) is continuous for all 
Lu and T as we hoped, but its T-derivative (essentially 
I{x,t)) is not. As T — > 0, the integrand of ([7]) becomes 
a spike of finite height and zero width. The integral past 
this spike is clearly zero (so the entropy would be zero 
as concluded above), but we run into trouble with the 
continuity condition. Were the limit T to be taken 
prior to integration, /(O, 0) would be zero instead and / 
would be continuous with respect to x but discontinuous 
with respect to t. 

Physically it does not matter whether /(O, 0) is 1, or 
something in between, since the contribution from this 
single point is zero in either case. Hence this one point 
should not matter. Formally we could state this by ex- 
cluding the point oj — Q from the uj integrals ([5]) and ([7]). 
Furthermore it is well known that the notion that every 
e{uj) on the form A + iBuj or A + iB/uj would violate 
Nernst's theorem is incorrect; on the contrary we argue 
that so long as A and B are temperature independent, 
none of these will. Using an analytical software such as 
Maple"""^ it is quick to check that integration of H{x,t) 



with respect to x followed by differentiation with respect 
to t gives the same result as when the order of the op- 
erations is reversed. While this argument is not rigorous 
it should convince the reader that the continuity issues 
at zero temperature and frequency can be avoided since 
this point is of no physical significance. 

A formally similar problem emerges when the permit- 
tivity is temperature dependent all the way down to zero 
temperature, as we will see, and in the latter case the sin- 
gularity at zero temperature and frequencies does appear 
to give a physical contribution and cannot be ignored. 



C. Uniform convergence 

An improper integral ([8]) is said to converge uniformly 
([lB| §4.42) if Ve > there exists a number oo > inde- 
pendent of T such that for all a, a' > ao, 



dxf{x,T) 



< e. 



Let us briefly analyse the behaviour of LnDq{u},k±) as 
\lu\ and k± approach infinity. The existence of the free 
energy integral ^ itself is well known, hence we need 
but check explicitly whether the integral ([7]) converges 
uniformly along different directions in the uj, ck± plane; 
clearly, if the double integral over co and k±_ converges 
uniformly, the oj integral ([7]) does so as well. 

As argued we consider propagating and evanescent 
contributions separately, in which case uniform conver- 
gence is straightforward to check. Reflection coefficients 
fall off rapidly as \uj\ — > oo (e.g. for the Drude model the 
real and imaginary parts of fall off as a;"'* and oj^^ re- 
spectively) and for > 1, kq is real and positive so the 
integrand furthermore decreases exponentially. The fac- 
tor exp(— 2Koa) is oscillatory for < 1, but the Dirich- 
let integral dxsinx/x is known to be uniformly con- 
vergent, and our integrand converges more quickly than 
this. It is easy to check that this also holds as \(3\ 
and \(3\ — s- oo. 

The splitting of ^ into propagating and evanescent 
parts may be done by integrating each part of the plane 
and taking the relevant limit to ck± —^ \uj\ in the end. 
Convergence problems are then avoided for the imaginary 
part of LnDq in ([7]); reflection coefficients fall off rapidly 
and further help is provided by the factor 



exp(a;/a'T) 



exp(— |w|/wt); |w|S>wt- 



The rate of convergence due to this factor depends on T, 
hence apparently cannot be used to demonstrate unifor- 
mity. We are interested only in low temperatures, how- 
ever, so by defining a finite upper temperature limit T 
above which the formalism is not valid, oq can be made 
T independent (dependent on T only) . The fact that the 
convergence of this factor alone is not uniform for infinite 
temperature is unproblematic, of course. 
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Thus we conclude that Nernst's theorem is satisfied 
for T independent e{uj) satisfying 1-3. The violation of 
Nernst's theorem in temperature dependent cases, as we 
shall see, can be understood as a direct consequence of 
violating the continuity criterion of Leibnitz' rule for im- 
proper integrals. 



III. TEMPERATURE DEPENDENT 
PERMITTIVITY 

In many models used in solid state physics, e is temper- 
ature dependent for all temperatures, and herein lies the 
source of much of the controversy over what is the correct 
theory of the Casimir force between plates of real mate- 
rials. The reader should note that the above theory only 
requires that permittivity be temperature independent 
for a finite temperature interval close to zero tempera- 
ture. Rather than rigorously generalising all of the above, 
suffice it here to discuss how the introduction of tem- 
perature dependent permittivity illuminates the entropy 
problems that emerge and hints at possible resolutions. 
In the following we will think physically in terms of pos- 
itive frequencies, bearing in mind that negative frequen- 
cies exert mathematically equivalent behaviour through 
the symmetry criterion [T] 

The models which have caused bother so far are the TE 
mode of ([5]) using the Drude model ([9]) when v{T) 
as T ^ (perfect lattice, no impurities) and the TM 
mode for dielectrics pU)) with ct(T) — >• 0. The two dif- 
ferent cases share many common traits, so analysing one 
of them in detail will suffice as illustration. Since the 
Drude model has already been treated in numerous ef- 
forts by both sides of the dispute (e.g. 0,1^1), we choose 
the dielectric for the below discussion. 

As should be clear by now, the troubles with entropy 
emerge for small frequencies at low temperatures. Let 
us from now on consider the only interesting frequency 
range in which a;^,7Ci; <C k\c^,ujQ, but making no as- 
sumptions about the relative magnitude of lj and ct/eo- 
Physically this corresponds to bringing uj and T close 
enough to the limit so that for all quantities which de- 
pend on their absolute values separately they may be 
replaced by zero, and only quantities which depend on 
their relative values, specifically the TM refiection coef- 
ficient, remains in question. In this case ttm with (HI 
inserted simplifies to 



w(eo 



(15) 



where v = LotQj a = x/ s where x and s are again suitably 
non-dimensionalised variables proportional to lo and a 
respectively. We have plotted the real and imaginary 
part of the squared refiection coefficients as is shown in 
figure[2]for illustration using trx, — 11-66 as reported for 
Si in 22]. 



We find that 3fJe{r^i^} is 



(eoo - l)V(eoo + 1)' 




0.5 



-0 . 5 




-0 . 5 




V 



0.5 



0.71 except for x <^ s (v k, Q) where it is unitary. Like- 



FIG. 2: Plots of the real and imaginary parts of r^M with 
(|10P for uo, y/LO"f <^ kj_c,iJo and ~ tj/eo as simplified in (|15|) 
with eoo = 11.66. Top: as function of independent x and s. 
Bottom as function of v = x/s. 



wise 3m{r^[^} for small s is approximately zero for the 
most part but increases to an extrcmum for small |t;| and 
thence decreases linearly through at w = 0, the same 
linear behaviour that led to 3m{(/)(ti;)} oc lo and the dis- 
continuity of figure [1] in the TM case before, which we 
argued was not essential. This time, however, there are 
additional discontinuities as s — > (equivalent to a — > 0). 
In particular SJmjr^jyj} (which contributes to the inte- 
grals) is everywhere except at a; = where it can take 
any value between its maximum and mininum (w ±0.079 
for €oo — 11.66). 

Now remember that the imaginary part of the squared 
refiection coefficient shown in figure [5] is to be multiplied 
with either the coth factor or its T derivative, both of 
which diverge as T/cu as uj 0. The result is an ex- 
ceedingly volatile behaviour of the integrands of ^ and 
([7]) near zero frequency and temperature, and the limit 
where both are zero can take any value between —oo and 
oo depending on the way the limit is taken! This con- 
trasts the bounded discontinuity shown in figure [T] in the 
temperature independent case. 

Furthermore, when e contains temperature dependent 
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parameters, entropy ([7]) will have an additional term 

TM 



duj 



e-^-°°cotli^ 

g=TE 1 



9(r 



dT 



(16) 



Additional entropy problems stem from this term. From 
(fT5|) one finds with a little algebra that 



d_ 

dT 



iv{eoo — 1) — 1 1 5(7 



(17) 



Assuming conductivity at low temperatures to behave as 
tJo exp(— Tq/T) with ctq a constant (see below), 



1 da 
adf 



y2 ■ 



When uj =/= this inverse quadratic temperature depen- 
dence is no problem since (jl7p varies as v^^ oc cr, so the 
term (fT6|) is zero by a good measure when T = 0. The 
limit Lu may be taken so that v has a finite value, 
however, in which case the derivative (jl7p diverges as 
T~^. This corresponds to r^j^ making a sudden jump 
from 1 to at T = for = 0. The term d{rl)/dT 
at T = is thus zero for all frequencies except uj — Q, 
where it is infinite. We can no longer argue that this 
one point does not contribute to the physical quantity S, 
and while the purported zero temperature entropy may 
be difficult to calculate in this way (it is straightforward 
to calculate it in the imaginary frequency framework), it 
seems likely that entropy obtains a finite value assuming 
the formalism may be taken at face value. 



A. Findings of the mathematical analysis 

An important realisation is thus that a violation of the 
third law of thermodynamics is predicted in the present 
framework when changes the power of its leading 
order term with respect to uj at exactly zero temperature. 
When t{uj) changes from oc to oc uj^^, a violation 
occurs in the TE mode; when the change is from oc lu^^ to 
oc oj^ , the TM mode gives the zero temperature entropy 
contribution. 

The bottom line is that when viewed in the frame- 
work of real frequencies, all apparent zero temperature 
entropy anomalies stem from divergencies of the Lif- 
shitz integrand in the double limit T ^ and w — > 0. 
The reader should note that in this author's understand- 
ing, this result does not contradict the findings of either 
[U, S H i, H, 01 or 0, i, i]. The latter of these refer- 
ences also notes a formal violation of Nernst's theorem 
due to rotation of permanent ionic dipole moments in 
the materials. In light of the above we may conclude 
that the latter effect is only problematic in media where 
the rotational degree of freedom of the ions vanishes with 
temperature in a manner so that its resonant frequency 
is finite at T > and zero at T = 0. This anomaly needs 
to be studied further in the future. 



Physical discussion of thermodynamical 
anomalies 



In this section we will undertake a brief physical dis- 
cussion of the mathematical results in the previous sec- 
tions, reviewing the temperature debate for metals and 
semiconductors in light of the above analysis. 

Models used when studying the physics of real sys- 
tems are founded on assumptions which we may cate- 
gorise as modelling idealisations and approximations in 
the description of the behaviour of these models, and 
there may be a need to distinguish between these in the 
present context. A relevant modelling idealisation, for 
example, is that a metal sample has infinitely large di- 
mensions and a perfect crystal lattice structure. Much of 
science is founded on such ideal models and corrections 
to them. A relevant approximation in this context is the 
use of simple dielectric functions such ^ to (fT^ which in 
particular assume that the media in question have local 
dielectric response (i.e. they only depend on frequency, 
not momentum k). Even for idealised systems, such ap- 
proximations typically have limitations. 

An ideal model which can in principle be realised 
(notwithstanding infeasibility in practice) cannot be al- 
lowed to violate the laws of nature, thermodynamics in 
particular. An approximation, on the other hand, will 
typically have a finite range of applicability, and cannot 
be expected to behave correctly outside this range. Given 
that the limits T = and oj = Q are in some ways ex- 
treme cases, it is especially important to investigate the 
latter point in relation to the purported problems with 
the third law of thermodynamics. 

Specifically, if an approximation which works well at 
room temperature does not hold for T = one cannot 
conclude from a formal violation of Nernst's theorem that 
it cannot be used within its applicability range. 

The much investigated temperature anomaly for met- 
als is a good example of the above, and we will review it 
briefly for illustration. For a perfect and infinitely large 
metal lattice, relaxation v{T) is due to electron-phonon 
interactions only and follows the Bloch-Griineisen for- 
mula (see appendix D of [3]), according to which v van- 
ishes as as temperature tends to zero leading to the 
above reported anomalies. It has been pointed out that 
no real metal sample is ever perfect [3| nor infinitely 
large [39}. so relaxation does not vanish in real systems. 
There is now consensus that for impure metals the Drude 
relation does not lead to thermodynamic inconsistencies. 

However, the theoretical problem is thereby only 
halfway solved, because as pointed out [2] the fulfilment 
of the laws of nature cannot hinge upon the presence 
of imperfections: the ideal system should accord with 
thermodynamics as well. The solution according to the 
authors of is to remove the relaxation from the Drude 
model and, more recently, introduce dissipation instead 
through a generalised plasma model [l6j , unfortunately 
at the expense of ignoring the manifest presence of re- 
laxatation at room temperature. Experiments seem to 
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confirm the predictions from such a procedure and rule 
out those impUed by the use of the Drude model © (e.g. 
psl p but the theory has not been universally accepted. 

While the ideal crystal lattice, when treated in all de- 
tail, should certainly be found to abide by Nernst's theo- 
rem, the approximation that its dielectric response is well 
described by the local formalism has been questioned for 
temperatures approaching T = 0. Svetovoy and Esquivel 
[26t and Sernelius [g^l conclude independently that at 
low temperatures nonlocal effects (the anomalous skin 
effect) dominate, and the local models are no longer re- 
liable. Their analyses accounting for spatially dispersive 
effects reveal that, within the approximations made in 
[2^ [131 Nernst's theorem is fulfilled independently of the 
presence of imperfections as it should. The spatial dis- 
persion approach was criticised [28j on several accounts 
with reference to a treatment by Barash and Ginzburg 
many years ago[29'| (see also (soj). The paradox remains 
that such a careful procedure (albeit not free of approx- 
imations) does not accord with available experimental 
data. Commendable efforts at a resolution include the 
recent exploits by Bimonte (3l| . 

While mathematically analoguous to the metal case, 
the temperature anomaly for semiconductors is physi- 
cally different. Here the problem is not related to ide- 
alised models (the conductivity of isolators truly does 
vanish at zero T), but approximations only. 

One can argue intuitively that the approximate model 
pn)) probably cannot be taken at face value when con- 
ductivity is very small since it implies that even an in- 
finitesimal conductivity should give rise to large thermal 
corrections in the Lifshitz formalism contrary to physi- 
cal intuition. For isolators, by definition, conductivity 
vanishes at T = 0, and for many such materials conduc- 
tivity even at room temperature is so small it would be 
expected to make for a minor perturbation only [i^]. If 
()10|1 is a poor approximation at low frequencies as a van- 
ishes, its violating thermodynamic laws in this case may 
not be too worrisome. 

A recent experiment measured the force between a 
substrate of the semiconductor silicon and a gold sphere. 
The semiconductor was excited into a metallic state by 
a pulsed laser and it was concluded that while the model 
was a good representation in the metallic state, the 
inclusion of the a term when the material was in the 
poorly conducting state was excluded at 95% confidence. 
This conclusion may not be surprising in light of the 
above argument which indicates that a Drude type per- 
mittivity model overestimates the effect of a small con- 
ductivity in the Lifshitz formalism. If so, it is likely that 
the experimental result might be explained without ref- 
erence to the Nernst theorem which concerns physics far 
removed from laboratory conditions. 

Attempts at a more careful description of the effect 



of a small density of free charges were recently rnade by 
Pitaevskii ^32,] and by Dalvit and Lamoreaux [33|, based 
on the effects of Debye-Hiickel screening from free charges 
in mean field theory. The resulting expressions do not fit 
experimental data well 3J| , and it is possible that a more 
detailed screening model is needed. 



IV. CONCLUSION 

While the Lifshitz formalism at real frequencies is 
much more complicated than the imaginary frequency 
equivalent normally considered, the consideration of 
quantities with direct physical interpretations sheds new 
light on the problem of non- vanishing Lifshitz entropy at 
zero temperature. We have argued that Nernst's heat 
theorem is not violated for any causal and continuous 
(except a,t bj — Q) e(w) which is independent of T near 
T = 0. This accords with the findings of Intravaia and 
Henkel[Tl| using a different approach. More generally, 
this holds for dielectric plates whose squared Fresnel TE 
and TM reflection coefficients are continuous for all ui and 
T and non-unitary except at cj = 0. It follows from this 
that entropy anomalies previously reported pertain to the 
persistence of the permittivity's temperature variance all 
the way to zero temperature and are consequences of di- 
vergencies in the Lifshitz formalism in the double limit 
w ^ and T ^ 0. 

When considering physical consistency in such limits as 
a means to distinguish between candidate theories, par- 
ticular care must be taken. We emphasise that approxi- 
mations can only be judged based on their performance 
within their domain of applicability. Specifically, approx- 
imations which are invalid at T = cannot be expected 
to be well behaved in this limit, hence cannot be rejected 
for causing violation Nernst's theorem which concerns 
zero temperature only. It is therefore important to ver- 
ify carefully that approximate physical models probed by 
invoking Nernst's theorem are valid in this case. We fi- 
nally argue that a recent experiment using an optically 
excited semiconductor can probably be explained with- 
out reference to the Nernst theorem by accounting for 
the presence of free charges more carefully than by the 
use of local Drude-type models. 
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